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ABSTRACT 

OPTIMIZATION OF HYBRID EVAPORATIVE COOLING AND 
AIRCONDITIONING SYSTEMS - AN ECONOMIC APPROACH 


A Thesis Submitted 

In Partial Fulfilment of the Requirements 
for the Degree of 


MASTER OF TECHNOLOGY 

by 

K. SEKAR 
to the 


Department of Mechanical Engineering 
Indian Institute of Technology Kanpur 

May, 1986 


A generalised computer programme has been developed 
for the cooling load calculations on hourly basis and the 
sol-air temperature variation for different orientation of 
walls* The effect of heat capacity and timelag factor of 
the structures on the cooling load estimation has also 
been studied. 

The single-stage vapour-compression airoonditioning 
system has been optimised to get the minimum total cost, 
for refrigerants R-12 and R-22. This is compared with a 



( xiii) 


hybrid system of evaporative cooling and airconditioning as 
to what extent, the latter is more economical than the 
former. The optimum results have been found out for the 
inside design conditions of 24°C, 65%; 27.5°C, 56% and 30°C, 
60% (dry bulb temperature and relative humidity) being 
obtained on the basis of Fanger's comfort equation [ 6] • 
Various combinations of the operating parameters are also 
studied. 
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CHAPTER 1 

introduction 


1 1 AIRCONDITIONING SYSTEM. 

Airconditioning is constantly absorbing and apply- 
ing the latest discoveries of science and engineering 
m an effort to ease man's environmental problems on 
earth and an space. The term/ airconditioning means the 
act of putting the air m the proper condition for the 
desired use. This may mean cleanliness/ the temperature 
and humidity of the air. Airconditioning for human 
comfort is becoming almost a prerequisite tor success m 
the merchandising field in all except the cooler regions 
of the country* The gain m customer's good-will is 
reflected by such an increase m business that the air- 
conditioning pays more than for itself. As any other 
engineering system/ the airconditioning system is controlled 
by various system parameters. To achieve maximum utilisation 
and effective implementation/ an air conditioning system 
has to be economical. In view of this/ the operating 
parameters of the system are optimised on the basis of 
cost as well as functional energy constraints whichever 
is contemplated. To carry out this objective/ there are 
alternatives available/ namely modification cf an existing 
system and design of a new system. The present work is 
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done on the basis of the former. 

1.7 SELECTION OF ENVIRONMENTAL ZONES rtND OUTDOOR 

TEMPERATURES 

Eventhough, there are many applications, the 
primary function of an airconditioner unit is considered 
to create temperatures artificially which will be comfor- 
table for human beings* Many investigators have done 
extensive work on temperature zones for human comfort. 

Fanger [ 6 ] elaborated the methods of getting comfort 
conditions and also suggested a comfort equation after 
carrying out extensive field studies. In India, Ramamoorthi 
[15] has analysed the equation for the Indian weather 
conditions after taking into account of several variables 
and suggested the comfort conditions, namely inside dry- 
bulb temperature, relative humidity and velocity of air* 

Singh and Prakash [14 ] presented a new ASHRAE comfort chart and 
recommended a region of design conditions for the comfort 
environment. The result of the field studies carried 
out by Malhotra [lO ] for the conducive atmosphere pertain- 
ing to Indian environment has also been taken into consi- 
deration. Based on the energy cost, Whitner [18] has got 
the higher comfort conditions. Several investigators [1 ] 
assured the daily to outdoor temperature variation as a 
sinusoidal function which results in Fourier series, in 
India, Ramamoorthi [15] suggested that the air conditioning 
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system can be substituted by evaporative cooling system, 
whenever feasible* 

The design of an aircond ltioning plant is mainly 
dependent on the cooling load, a function of indoor and 
outdoor temperatures and internal energy release. The 
indoor temperatures are chosen on the basis of comfort 
industrial requirement. The size of the building and 
various capacities of the equipments are also the primary 
factors influencing the cooling load calculations. So, 
great care should be taken in proper selection of the 
equipments to get the total cost as minimum. It has been 
a common practice to use air-cooled condensing unit for 
both residential and commercial applications. This 
dispenses with water piping, cooling water tower and 
water treatment, thereby eliminating their costs. But the 
head pressure »f the air-cooled system is usually 10% to 
20% higher than that of the water cooled system. Hence 
the air-cooled system is not used in case of large cooling 
requirement. For smaller units /the cost and conveniences 
dominate over other factors* Hence, the air-cooled 
system is considered in the present analysis* 

1.2.V ANALYSIS OF COOLING LOAD CALCULATIONS 

There are number of parameters affecting cooling 
load calculations/ often difficult to define precisely, and 
always intricately interrelated. Most of these parameters 
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can be assembled under the classification of space heat 
gain. This is necessary because different fundamental 
principles and equations are used to calculate different 
modes of energy transfer. The major parts of the space 
heat gain are (a) heat transfer through the structure of 
the building^ (b) heat generated within space by occupants/ 
lights and appliances and (c) ventilation and infiltration 
of outdoor air • For the evaluation of structural heat 
gain/ there are different techniques available, namely 
(a) sol-air temperature method/ (b) cooling load temper- 
ature difference method, (c) finite difference method and 
(d) finite element method. Among these, sol-air temper- 
ature technique is mostly preferred due to its simple 
algorithm. The cooling load temperature difference 
method [ 3] is implemented by using the transfer function 
method to compute the one -dimensional transient heat 
flow through the building structure. Kadambi and 
Hutchinson [9 ]have described the finite difference method 
applied for the heat transfer through the walls. Huebner [7] 
has elaborated the finite element method xn the applica- 
tion of heat flow through structures. The thermal capacity 
and thickness of the structure are interrelated and 
presented m the graphical form [ 2 ] . The relation between 
the infiltration load and the volume of the conditioned 
space is explained -in [3,4/12 ] * Woods et al [19] have 
modified the ventilation standards by giving emphasis on 
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energy saving The heat generation within the conditioned 
space has been elaborately given m [12 ] taking into 
account of metabolic ratings. The procedure for hourly 
cooling load calculation is elaborated m [3 ] . 

1.2.2 METHOD OF OPTIMISATION 

In general /optimisation can be defined as the act 
of obtaining the best result under the given circumstances. 
In design, construction and maintenance of any engineering 
system, engineers have to take many technological decisions. 
The ultimate goal of all such decisions is to either 
minimize the effort required or maximise the desired 
benefit . Since the required effort or desirable benefit in 
any practical situation can be expressed as a function of 
certain design variables, optimisation, in particular, 
may be considered as the process of rinding the conditions 
that results in the maximum or minimum value of the 
function, known as objective function. It is very rare 
that a practical design problem to be unconstrained. 
Considering this fact, if the expression for the objective 
function and the constraints are fairly simple in terms 
of design variables, the classical methods of optimisation 
technique can be used to solve the problem. On the other 
hand, if the problem involves the objective function and 
the constraints which are not stated as explicit functions 
of design variables or which are too complicated to 
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manipulate, we should go m for some specialised numerical 
methods of optimisation. The basic philosophy of these 
methods to produce a sequence of improved approximations 
to the optimum value. In the present work, one of such 
methods called as 'Penalty Function Method' [16] is 
adopted to obtain the optimum design variables of the 
airconditioning system for the minimum total cost. The 
application of this iterative procedure is made possible 
due to the availability of the computer. These methods 
are clearly elaborated in [16] • 

1 . 3 PRESENT STUDY 

In the present study, a generalised computer 
programme has been developed for the accurate evaluation 
of the cooling load considering daily outdoor temperature 
variation based on actual as well as approximate. In 
this regard, the effect of timelag factor and heat 
capacity of the structures on the structural lo^d has also 
been found out. The cooling loads have been calculated 
on the basis of actual hourly daily temperature variations 
and the sinusoidal temperature variation based on maximum 
and minimum temperatures. The multiplier has been developed 
to facilitate the utilisation of the suggested sinusoidal 
temperature variation; because these data are available 
for most of Indian cities [11]* 
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For the iterative optimisation technique/ a 
versatile computer programme has been developed which 
takes care of various factors involved in the total cost 
of a vapour compression system. Even though there are 
large number of operating variables present m an air 
conditioning system/ the most important factors have been 
selected as design variables/ giving emphasis on the 
feasibility of control. For these variables/ the constraints 
are formulated on the basis of practical consideration* 

A comparative study has also been carried out between the 
airconditioning systems operating with and without evaporative 
cooling/ for the duration/ April to June/ in terms of total 
cost. This aspect is studied for both R-12 and R— 22 
refrigerants. 
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CHAPTER 2 

COOL IN 3 LOAD ESTIMATION 


2.1 INTRODUCTION 

The design of an airconditioning plant for a 
building depends on the use of the building. The cooling 
load calculation needs the size of the building/ its 
structural details and orientation. In addition.* the 
following cfuantit le s are considered. 

l) Design conditions* Indoor and outdoor conditions* 
n) Instantaneous heat load* Sensible and latent. 

For the given inside design condition/ the cooling loads 
are calculated as follows. 

Sensible cooling load It comprises the heat 
transfer through structure including solar heat gain. 

The sensible heat part of the infiltration and ventilation 
is calculated from the known amount of the air and subsequent 
difference between enthalpies of the inside and outside 
air. The sensible heat part of the heat release from 
the occupants is calculated from the standard data 
available giving heat release from human beings* ratings 
of fans* lights* etc. 

The latent heat part of the cooling load has been 
calculated from the Infiltration and ventilation air* 
latent heat release from the occupants Inside the room. 



9 


2.2 SELECTION OF INDOOR CONDITIONS! 

The variation m the weather# building occupancy# 
and other factors, affecting the load# necessitate 
carefully coordinated controls to regulate simultaneously 
the components and equipments to maintain the desired inside 
conditions. The design indoor conditions play a key role 
m the design of any comfort air conditioning system. The 
field studies [lo]# reveal that m India# for hot humid 
climates# the comfort zone is between 22°C and 24,5°C# 
effective temperatures. The analysis [15] carried out on 
the basis of Fanger's comfort equation recommends the 
comfort conditions for Indian people as T db = 24°C# 

RH = 0.65# V = 0.8 m/rnin, T db = 27.5°C# RH = 0.56# 

V = 1.1 m/s, and T^ b p= '30°C# RH = 0.60# V = 0.7 m/s. 
Since# these results lie m the comfort zones recommended 
in [5#6,i43#in the present analysis the same combinations 
are taken into consideration. 

2.3 HOURLY COOLING LOAD. 

Most of the factors of the cooling load vary in. 
magnitude over a wide range during a twenty four hour 
period. As the cyclic changes in the load components are 
not usually in phase with each other# the analysis is 
required to establish the resultant maximum cooling load 
for a building. 
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2.4 HOURLY VARIATION OP OUTDOOR CONDITIONS 

Hourly variation in the outdoor conditions is 
available for the limited cities. But the cooling load 
depends greatly on the ambient temperature variations 
as compared to that obtained on the basis of design 
conditions. So a mathematic model [1] has been adopted 
which gives a close approximation to an average value 
equal to the mean daily temperature plus an oscillating 
component with an amplitude equal to one half the mean 
daily range. For most practical applications/ this can 
further be approximated by a constant with a superimposed 
sine wave having a 24 -hour period. So, the final form 
of the- -expre-ssiojo. comes out to be a harmonic series of 
first order 


T = A + B cos ( 15 t - C) 
o 


( 2 . 1 ) 


where T 

o 

t 

A/B,C 


outside air temperature at any time, °C 
time m hours 

the constants being obtained by applying 
the boundary conditions/ 


dT 

o 

dt 


= 0 when 


(*T = T 

{ 

\T^ = 

. .s o. 


( 2 . 2 ) 


min 


T / occurs at (Sunrise-1)'' ^hrs? T” 1 *** occurs 12 hrs thereafter* 
min max 

So, the above model is solved for the outdoor air temperature 
variation by knowing the minimum and maximum temperatures 
with their corresponding timings/ in a day. The result of 
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this equation has also been verified with actual variation 
of the temperature and the comparative study has been 
done in a graphical form. Fig. 4,1, for Kanpur. For this 
purpose, the data for the actual variation of the temper- 
ature and relative humidity for the period 1982, 1983 
and 1985 have been procured [8] and given in tabular form 
in Appendix-A. 

2.5 SOLAR RADIATION. 

The solar radiation is the sole factor for the 
summer air conditioning system. The total incident solar 
radiation is given by 

It = I DN Cos© + I d , k W/m 2 (2.3) 

where the product of Ij^, direct normal radiation and the 
cosine of the incident angle ©, stands for the irradiation 
on a surface from the sun and 1^ refers to the diffused 
radiation. Simplified general relations for 1^ and 1^ 
for a clear sky are given approximately by, 

2 

Id N = A 1 /exp (B-j/SmP), k W/m 

J d “ <V • F ss ' 15 W/m2 

where, 

A^ = apparent solar irradiation at an air 

2 

mass =0, k W/m 


2.4(a) 

2.4(b) 
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= atmospheric extinction coefficient , 

C x = diffuse radiation factor, 

F gs = the angle factor between the surface and 
the sky, 

= 0,5 for vertical surfaces, 

=1,0 for horizontal surfaces, 

= (1 + Cos$J^)/2 for any other inclined surfaces 
and P = the altitude angle. 

The values of A^, and C 1 of [3 ] are given in the 
Appendix -B. 

2,5.1 ESTIMATION OF SOL-AIR TEMPERATURE 

For heat transmission calculations, it is convenient 
to combine the effects of outside air temperature and 
solar radiation clustered into a single fictitious quantity 
called sol-air temperature (T gQ ^). The rate of heat 
transfer, q, from outside to inside surfaces of a sunlit 
structure may be written as, 

<3 = h o (T o - T w ( o> + *t • < 2 * * - 5 > 

where, 

h - outside air convective heat transfer 
o 

2 

coefficient, k W/m - °C 

T q = outside air temperature, °C 

T * temperature of the outside surface, °C 
w,o 
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a^ = absorptivity of the structure/ 

/ 2 

and I^_ = total intensity 01 solar radiation/ Jc W/m • 

The total intensity of radiation is calculated by the 
primary angles namely/ the latitude (1)/ hour angle (a j^) 
and the sun's declination/ d. Besides this/ there are 
several other angles as the sun' s zenith angle ¥ , 
altitude angle (3 / and azimuth angle Y « In Appendix-C, 
the equations for the estimation of fchese angles are 
given* 

The equation (2.5) can be modified as 


g = h (T . - 
H O sol 

T ) 

W/O 

( 2,6) 

where , 

T i = T + 

sol o 

a «r I /h 
t t. o 


The outside air convection 

heat transfer coefficient. 


h^is calculated as follows 

h = 14. 27 8+3. 5 59 .V 

O w 

for very smooth surfaces. 

2.7(a) 

= 18.423+3. 81. V w 

ior smooth surfaces. 

2.7(b) 

= 26.55+5.0663 .V 

Vjr 

for rough surfaces. 

2.7(c) 

= 28.64+6.364 .V w 

for very rough surfaces. 

2.7(d) 


where V is the design velocity of the outside air m 
w 

Km/h. In the present study, the variations in the sol- 
air temperature for a day/ on different orientation of the 
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surfaces presented in a graphical form. Fig. 4. 2. 


2.6 HEAT TRANSFER THROUGH STRUCTURES. 

A fraction of the total heat energy is absorbed 
by the structural material and there is also a time lag 
for the heat transfer from the outside surface to the 
inside surface of the structure. The thermal capacity 
of the walls and ceilings/ taken as/ 


c th - m. c p = P . A . c p (2.8) 

where rn, c^/ P , A and x are massy specific heat / density , 
cross-sectional area and thickness of the structure/ 
respectively. 

Taking into account of timelag and decrement factor 
the actual heat transfer through the structure at any 
time l t* is given by 


♦ 



5 - 

S [Uj 

j=l 



(T om " V + Yj*V T o(t-t) 



C 2.9) 


where T , A- and T are the mean sol-air temperature/ 
om' 

decrement factor and timelag factor/ respectively. In the 
above expression, u, A, and T_^ t) re ^ er 
heat transfer coefficient, cross sectional area and 
temperatlww at time (t-T), respectively. Here summation 
symbol stands for four walls and the ceiling- 




K, =0 00072 kW/m°C 
K 2 =0 00130 kW/m®C 


hi |T ( 


Fig 2 1 (a) Structural details of the wall 


K 3 ==0 00072 kW/m°C 
K 4 *0 00173 kW/m # C 

Kj 

K 3 

Fig 21(b) Structural details of the ceiling. 
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Referring to the Fig. 2.1, the heat conduction 
through structure can be written as/ 

(T, 


w ,o 


^st h o (T o T w / o ) i =3 


- T ) 


U (»/K ._) 

1=1 


( 2 . 10 ) 


The above equation can also be written asy 


l st 


= U 


(T o * 


T i) 


( 2 . 11 ) 


with 


U = 


1 



1=3 

2 

i=l 



T q = outside air temperature y °C 
T = indoor temperature/ °C 


x and K, 
i i 


h 


i 


u 


thickness and thermal conductivity of the 
'i‘ th layer of the structural material 
inside convective heat transfer coefficient/ 
kJ/m 2 - h - °C 

overall heat transfer coefficient/ kJ/m -h— °C, 


In the present study, the values of K 1 a are taken from [3]« 
The inside convective heat transfer coefficient is found to 
be function of temperature difference (At = T ^ - T^). 
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For walls? 

h A =4.25 (A T) l/4 , kJ/h - m 2 - °C 2.12(a) 

For ceiling 

h x = 3.48 (At) 1>/4 , kJ/h - m 2 - °C 2.12(b) 

The Newton-Raphson iterative procedure has been used to 
get the values o£ ' h^ 1 • Equation 2.9 is used to calculate 
the heat transfer through the structure after considering 
time lag and decrement factors [2] • The effect of these 
factors on the cooling loads has been discussed, m 
Chapter 4. 

2.7 heat transmission through glass windows 

The incident radiation may be reflected/ partly 
absorbed and the remainder transmitted through the glass 
materials. Thus, m general, it may be written as 

r + a_ + t = 1 (2.13) 

g g r 

wnere r , a and t are the reflectivity, absorptivity 
9 <3 r 

and transmittivity, respectively. The equations for 
calculating these thermal properties for the glass windows 
are given m Appendix -D. 

From the basic concept of heat transfer, the rate 
of heat gain through the glass is given by. 
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g glass F s* t rD* I D +t rd* ;r c 3+t 


_.I D +h CT - T ) 
rR R i,g g,i i 


( 2.14) 


where P g is the sunlit fraction of the window surface, 

^rd 7 an< ^ t r R are the transmittivities of the direct, 
diffused and reflected radiations, respectively. 1^,1^, 
and I^ are the intensities of direct, diffused and 
reflected radiation, T and T. are the temperatures of 
inside surface of the glass and indoor conditions, 
respectively. 

In reality, a glass window is not fully exposed 
to direct sunlight, but partially shaded. Referring to 
Fig* 2.2(a), it can be written as . 


X = b 

Y = b 


i tan a 

tan 6 

Xr 


2.15(a) 

2.15(b) 


where the profile angle 6 is related to the sun's altitude 

XT 

angle p and wall solar azimuth angle a by. 


tan 6 p 


tanP 
Cos a 


( 2.16) 


Finally, the value of sunlit fraction»F t is given by, 

s 

F g - l-P^tan^ - P 2# tana+ p^. 9^, tan tan a (2.17) 

where , 

P^ = b/a and P^ - b/c. 
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Referring to Fig. 2. 2(b), the energy balance for the glass 
sheet can be written as. 


s * * a gD + ^d * a gd 


+ I„.a _ = h (T -T ) + 
R gR i,g g,i i' 


h (T -T ) 
o,g g,o o' 


( 2.18) 


where h. 


an ^ ^o,g are ^e sur ^ ace coefficients given by 


h o,g * 3,96xlO" 3 + 9,89 2xlO~ 4 .V o , kJ/m 2 -°C -S 2.19(a) 

h 1/g “ 0.237+( V^/Height) 0 * 5 , kJ/m 2 -°C-k 2.19(b) 


where V Q is the outside air velocity in Km/h and V is the 
indoor air velocity in m/h. 

Combining the equations (2.14) and (2.18), the heat transfer 
through windows and doors including solar heat gain is 
calculated as 


N 


°9 “ ^ =1 + fc rd * X 6 + t rR I R ) + 


! V „ „ (T . T) , 

1 TT _/h. ) + U \ T - T ± ) j 

o,g i,g 


( 2 . 20 ) 


where. 


U = 


(> — + 
i,g o,g 


N 


number of windows and doors 
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The details of windows and doors are given m Appendix-E. 
The total heat transfer through the structure is given 
by the sum of equations (2 9) and (2.20) as: 


'st 


= Q t + 


( 2 . 21 ) 


2.8 INFILTRATION LOAD 

Infiltration is the leakage of outdoor air into 
a building through cracks and opening caused by pressure 
difference across the boundary surfaces. The net exchange 
of air may lead to both heat and moisture gain for the 
space The volume of infiltration is related in terms of 
the room volume and taoular values are available m [3]. 
The expression for the infiltration load is 


V 


'inf il 


room 


air 


51 N ACH (h a - h x,a )/24 ' k W 


( 2 . 22 ) 


where V is the volume of the room/ v . , is the 

room 

specific volume of air at outside temperature/ 13 

the number of air changes/day, h a and h^, are the enthalpies 
of outside and inside air, respectively. 


2 .9 VENTILATION LOAD. 

The introduction of outside air into the conditioned 
space is an important factor, which is necessary to maintain 

and odour levels as per standard practice. Ventilation 


oxygen 
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codes and standards recoituxiend a minimum rate of 30 cfm/ 
person. But the recent analysis shows that it is more 
efficient at 5 cfm/person [19] on the basis of economic 
consideration. Since the present problem deals with 
economic aspect, this value has been adopted for the 
ventilation load calculation, which is given by 

°vent±l “ * ( ventil/person) * ( h. - h > 

air 

( 2.23) 


where v_ is the specific volume of air and h , h. „ are 
air a 

the entnalpies of outside and inside air. 

2.10 MISCELLANEOUS LOAD: 

Apart from transmission, there is a considerable 
amount of heat energy generated inside the conditioned 
space from the occupants, lights, electrical equipments 
and appliances. The thermal energy release frcm human 
body at different degree of activities is taken from [l^'^L 
In the present analysis, the degree of activity is taken 
as the average cf all activities and the corresponding 
heat release has been calculated. The instant rate of 
heat gain from lights and other electrical appliances are 
calculated as. 
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Q e iec = ( Number o£ electrical appliances x rating of 

the appliances) + (Number of lights x ratings) x 
cooling load factor 


( 2.24) 


In the present work/ the cooling load factor is taken as 
0.88 [15] . So the miscellaneous load is calculated as 


mi sc 


= Q 


elec 


occupants 


( 2.25) 


The total cooling load is calculated from 


^total ^st + ^infil ^Ventil + ^misc 


( 2.26) 


Using this procedure/ m the present analysis, hourly 
cooling load calculations have been carried out for 
a building model. Fig. 2.3, for Kanpur. Fcr this estimation, 
a generalized computer programme has been developed with 
great care and is given m Appendix-G. 
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CHAPTER 3 

PROBLEM FORMULATION AND OPTIMIZATION 


3.1 INTRODUCTION 

The problem formulation for optimization is done 
on the 1 »asis of total cooling load based on the design 
condition for the equipment selection. The total cost 
calculated on the basis of fixed and running costs for 
comfort airconditioning constitutes the objective function 
m terms of the operating variables having a set of 
constraints. The optimization problem can be stated as 

which maximizes or minimizes f(X) subject to 

■+ 

the constraints g m (X)^ 0, m = 1/2/. .ry where X is an 
n -dimensional vector called as the design vector / £(X) 
is the objective function and g m (X) are constraints • 



3.2 DESIGN VARIABLES. 


Any engineering system is described by a set of 
quantities some of which are viewed as variables during 


the design price ss. In general , certain quantities are 
usually fixed at the outset, known as preassigned parameters 
and all the other quantities are treated as variables in 


the design process and collectively represented by the 

"* . Even though there are many variables which control 


vector X 
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COMPRESSOR 



(b) 


Fig 3 1 (a) Single-stage vapour-compression system 
(b) P-h diagram for single-stage vapour - 
compression cycle 
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the airconditioning performance, inclusion of all these 
variables will lead to enormous computational effort. 
Here, an optimal prediction effort has been made to 
select the variables judiciously. 


3 . 3 SYSTEM DESCRIPTION 

In the present study, a single-stage vapour 
compression refrigeration system is taken for the 
analysis having refrigerants R-12 and R-22. The cycle 
is shown schematically m Pig. 3 1# To simulate the cycle 
close to represent actual processes, a pressure drop m 
the compressor valves, heat rejection in the compression 
discharge line and the degrees of subcooling of condensate 
and superheating of evaporator vapour are considered. 

The pressure drops in the condenser, evaporator and 
piping are ignored. In the present analysis, air-cooled 
condenser has been incorporated m the system. 


In the design of the cycle, the pressure drop at 
the suction line is assumed as 0.2 bar [17]* The 
increase m pressure at the delivery of the conypressor 
is taken as 0.4 bar [ 17] . In the present optimization 
analysis, four variables are considered. They are. 


xu> 

X( 2) 


Degrees of rise xn the condensing temperature. 



Degrees of fall m evaporator temperature. 
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X( 3) Degree of superheating/ T , # C 

s 

X( 4 ) Degree of subcooling/ T c /°C. 


3.4 ESTIMATION 01 OBJECTIVE FUNCTION. 

In the present study/ the objective function is 
taken as the total cost for the comfort airconditioning 
system# which combines the initial and running costs of 
the system. If C I and C R are the initial and running 
costs per year respectively/ then the total cost/year 
(Qp) is given by 


Op = C I + C R , Rs/year 


(3.1) 


where 


and 


C I = 


C comp + G evap + G cond + Slower 


C = Cost of the compressor# Rs/year 
comp 

C = Cost of the evaporator# Rs/year 
evap 

C , = Cost of the condenser# Rs/year 
cond 


Si 


ower 


=• Cost of the blower 
= C 1 + c 2 + C 3 + C 4 


Rs/year 


3.1(a) 


3.1(b) 


Cp = Cost of the compressor power# Rs/year 

q = cost of the blower power# 


= cost of the fan power# 
= Maintenance cost# 


Rs/year 

Rs/year 

Rs/year, 


The 'Present Worth Method' is used for the analysis of 
the economic model. 
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3.5 OPTIMIZATION TECHNIQUE 

The presence of constraints in a nonlinear progra- 
mming problem creates more problems while finding the 
minimum. Several situations can be identified depending 
on the effect of constraints on the objective function* 

The simplest situation is when the constraints do not 
have any influence on the minimum point. From the 
application point of view on the real problems, there is 
a possibility of having two or more local minima for 
the minimisation of the problem. This situation is entirely 
due to the nature of the objective function contours. By 
taking into account of all these possibilities, in the 
present analysis, one of the indirect methods of non-lmear 
programming constrained optimization technique, known as 
'Interior Penalty Function’ [ 16 J is used, 

3.5.1 THE PENALTY FUNCTION APPROACH S 

Penalty Function Methods transform the basic 
optimization problem into alternative formulations such 
that numerical solutions are sought by solving a sequence 
of unconstrained minimization problem. For analysis 
purposes, a basic optimization problem is taken such that 
X minimizes the objective function, f(X) subject to the 
constraints g m (X) < 0 ,m = l,2,., *ry This problem 
is converted into a unconstrained minimization problem bys 
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0 k = 0(X, r^) = f(X) + r fe S 1 G. [g(X) ] (3,2) 

j=l 3 

■4 

where Gj is the function the constraints/ g( X) / and r^ 
is the positive constant known as the 'Penalty Parameter'* 
The function/ 0(X , r^) / is optimised for a decreasing 
sequence values of the penalty parameter/ r^. Evidently, 
if the 'r^' is vety small/ the optimization of 0(X, r^) 
amounts to optimization of the objective function/ f(X). 

3,5.2 INTERIOR PENALTY FUNCTION METHOD: 

In this method, the function Gj[g(x)3 of the 

equation (3.2)/ is substituted by - ■ — • The 

<x) 

unconstrained minima of 0(x y r^) / all lie in the feasible 
region and converge to the solution of f(X) as r^ is 
varied in the decreasing sequence. So the function 

4 

0(X, r^) is defined as/ 

n l 

0( X, r.) = f(X) - r. 2 ks — (3.3) 

k k m=l g m (X) 

Here/ the values of function 0 will always be greater 

4 

than f/ since g m (X> is negative for all feasible points/ 

X. Since the unconstrained minimization of 0(X, r^) is 
to be achieved by changing the sequence of penalty 
parameter/ r k , in the decreasing order, its initial value 
should be chosen such that quick convergence is obtained* 
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Thus, for any feasible starting point/ the value of 
r^ can be taken as/ 


0.1 f 

(X x ) 

f (X 1 ) 

v — ^ “ 

. & 

1 

s r- <* JL 

1 £1 1 


- ^ 

m=l 

9m <V 

m=l g m ( Xl ) 


(3*4) 


Once/ the initial value of r^ is chosen/ the subsequent 
values of have to chosen such that 

r k+i < r k 3 - 5(a) 

In this regard/ for convenience purposes/ it is taken as/ 

r k+1 « D r fe/ where 0<D<1 3.5(b) 

where in the present analysis/ D has been taken to be 0.1, 

In the present analysis/ the unconstrained 
minimization of the penalty function 0(X, r^) is done by 
'Davidon-F letcher-Powe 1 Is Method' [16]. This is the best 
general purpose optimization technique making use of 
the derivatives that is currently available. The algorithm 
for this iterative technique is as follows: 

(a) Start with an initial feasible point x^ and 
a nxn positive definite symmetric matrix /[HJ which 
is taken as [i] • 
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( b) Compute the gradient of the function/ 7 0 

at the point X 1 and set = - [H^J . V0^ . Normal sing 

s. 

this, we obtam> S, T = where 'n' is 

N1 i 1/2 

the number of variables. 

In the present formulation, central difference 

scheme has been adopted for gradient evaluation. 

★ 

(c) Find the optimal step length T m the direction 

*“► -V * 

of S and set X t , = (X + T S. T ) . In the present 
Ni 1+1 1 i Ni 

★ 

study, the determination of T is done by 'Cubic Inter- 
polation Method' [16]. 

(d) Test the new point X , for optimality using 

Cv h + i - 

the condition rf fl — E 1 ' w ^ ere e ± ls ver Y 

small. In the present case, the termination criteria 
—3 

is £ -j. < 10 . If this criteria is not satisfied [H] 

is updated as follows 




iamhj+p*,.:! 


where. 


M. 


T l • S HJ. ‘ S H1 

<4 • v 


N, 


-< H A> AA > 


v J* A+i> - 7 >*V 
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[ 


START WITH AN INTERIOR FEASIBLE POINT tT 

, -i 

|SET i>1 , K*1 I 

I 


CALCULATE THE VALUES OF OBJ- 
ECTIVE FUNCTION AND CONS - 
TRAINTS GG CM) 

.*U/GG(M) 
/06(M) 


CALCULATE 0(X, R k ) » f (7 ) -R k IL"1 
HERE, R» IS TAKEN AS -CXOBJ FUN 
WHERE, 0<C<1 
L 1 


FIND #(VX|) «V*(X',) SET Si * -[H,]x 
WHEN f«1, [H]»[I] 


|SET 7,+< » % »(T| )sT WHEN i » T, T ■ 1 1 

Find ~| 



YES 




jFlND Tf by CUBI^ INTERPOLATION METHOD 


[CALCULATE [M,]«T,*(S,KS^/(S| T )(a| ) 
[N, 3 - -(H,Q|)(H|Q|) T /H| q, t 

I 


|set [Hki3-IHiHm,]*[Ni3 1 
I Hr 



SET ftiM-tCHRK) 


^OPTIMUM *^•♦1 

WHERE 0< C<1 0 


t optimum ** ft ♦! 


H SET R K « Bk*1 1 


Fia 3 2 Flow chart for the interior penalty function method 
using Davidon- Fletcher - Powell's unconstrained 
optimisation technique 





ENTER WITH X, SAND T 0 

HERE T ft «1 0 

— - ♦ 

SET t A » f (0) , f(0). li AND A«0 




Fig 3,3 Flow chart for cubic interpolation 
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(e) Set the iteration number as i = i+l and 
repeat the procedure from the step 'b'/ until convergence 
is achieved. Flow charts for 'Interior Penalty Function 
Method' and 'Cubic Interpolation' are schematically presentf^ 
m Figs. 3.2 and 3.3, respectively. 

3.6 ESTIMATION OF LOADS ON EACH COMPONENT OF THE SYSTEM. 

The components of a refrigeration system have to 
be selected on the basis of the design conditions* The 
power requirements of the apparatus are evaluated for 
the design conditions and a factor of safety is multiplied 
to ensure the capacities of various components from 
the safety point of view. The condensing and evaporator 
temperatures of the vapour -compression systems are considered 
to be the sole factors which influence effectively the 
design of the airconditioning system. The evaporator 
temperature is selected on the basis of the indoor conditions 
whereas the condensing temperature is selected on the 
basis of design temperature. 

So the condenser temperature is taken as, 

T. = T + A T + X( 1 ) (3.6) 

h o * 

where/ 

T s= the dtesxgn outside air temperature/ C 
At » the approach temperature of the condenser, °C 

3 . 

X(l) = degrees of rise ±n condensing temperature, 
taken as a variable, °C* 
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Fig 3,4 Experimental set-up of an airconditioner. 
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The evaporator tenperature is chosen on the basis of 
inside room temperature. So it is taken as. 

T, = T - AT - X( 2) (3.7) 

1 i a 

where, 

T = the inside room temperature, °C 
At = the approach temperature, °C 

O 

X( 2) = the variation in tfte lower temperature 
which is taken as a variable, °C. 

In the present analysis, the values of the approach 
temperature have been^as 6°C, 10°C and 14 °C* To justify 
these values / experiment has been done on an air conditioner* 
The pressure gauges were fitted m the suction and discharge 
line. Fig. 3.4. The airconditioner has been operated and 
temperatures of air leaving the condenser and evaporator 
were taken. The air was supplied to both at about 30 C. 

The average aPi roach temperatures for the condenser and 
evaporator sides were found to be about 10 C. 

The value of T ^ is also checked up m the basis 

of the apparatus dew point temperature, T a(3p . Because 

the lower temperature, T^, should alvays be less than 

m <rhe T . is calculated in terrtls of effective sensible 

adp adp 

heat factor, ESHF, which is given by 
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ESHJ 


1 026 (T - T , ) 
x adp 


1 026 < Tl - T adp ) + 2410 (W a - W adp ) 


(3.8) 


where , 


ESHP 


ERSH 

ERLH 


W 


l 


W 


adp 


ERSH 

(ERSH + ERIn) 

RSH + (BPP) . OASH 
RLH + (BPP) . OALH 

Specific humidity inside the room, kg/kg of dry aii^ 

Specific humidity at apparatus dew point 
temperature, kg/kg of dry air* 


RSH and RLH represent the room sensible and latent heats 
respectively. OASH and OALH denote the outside air sensible 
and latent heats, respectively. The effect of the bypass 
factor (Bpf) on the performance of the cooling apparatus 
is discussed m [4"] and 0.15 is the value recommended for 
row coil. By using Newton-Raphson iterative procedure, 
the value of T ad p, is obtained, for the known values of 
ESHF. 


Referring to the Fig. 3.1, after taking the degree 
of superheat X(3), the temperature at the state point 'l* 
is given by. 


T x = T x + X(3) - A T x ( 3.9) 

where At^ is the decrease m the lower temperature befo-e 
compr e s sion , 0 C. 
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Similarly, T 2 = T h + ^ vT h 


(3.10) 


where A is the increase in the condensing temperature 
after compression, °C, 

The enthalpy values, h p f the refrigerants, R-12 and R-22, 
have been computed from the functional relations of the 
temperature of the state. These relations are given in 
the Appendix~F, 

From the Fig. 3.1, - X(4) (3.11) 


where X(4) is the degree of subcooling. 


The mass flow rate of the refrigerants, 
from 


m ref 


Q c x SF 
- h 4 ) 


'c. 


m 


ref' 


is found 


( 3.12) 


where Q c is the cooling load as given in Sec .2.10 and 
SF is the factor of safety taken as 1.25 [4]. 

The compressor power is given by. 


m re£ 1 h 2 - h l> 


comp 


T) 


, kw 


m 


( 3.13) 


where Ti is the mechanical efficiency of the compressor 

Tti 

taken as 0*85 [12} • 

For the air— coo led condenser/ the load is 


Q . e ni - (h 0 - h.) / kW 

w cond ref 2 3 

having h 2 = + ( A h) lsentr ppx</ n c * 


( 3.14) 
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The equation (3.13) can also oe equated to the heat gain 
by the cooling air as 


Q 


cond 


m 

air 


* C P 


air 


(T 


o,c 


T 

air ,1 


) ,kW 


3 . 15( a) 


where , 

m a ^ r mass flow rate of cooling air, kg/s 

c specific heat of air, kJ/kg.°C 

p air 

T the outlet temperature of the air,°C 

0,0 

T a ^ r x the inlet temperature of the air,°C. 


In the a cove equation, the temperature differential is also 
equal to the variation m the condensing temperature (T^) 
denoted by X(l). So the equation 3.15(a) can be rewritten 
as. 


, = m x c x X(l), kW 

cond air p air 


3.15(d) 


By combining equations 3.14 and 3.15b, 


m 


m re£ (h 2 " h 3 } 


air 


X( 1) 


/ kg/s 


(3.16) 


air 


Prom equation (3.16), the volume of the air handled by 
the fan for the condenser is calculated as. 


Vol 


air 


m, 


air 


x 60 


air 


irfVmin 


(3.17) 


3 

The fan power is correlated with its capacity (m /min) 
in the form of a quadratic equation. For the purpose, a 
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/IBHUMIDIFICATION air 



F - FAN 

S - WATER SUMP 



(a) 


A - OUTDOOR DESIGN CONDITION 



(b) 


Fig 3.5 (a) Evaporative cooling arrangement in rooirg 
(b) Evaporative cooling process. 
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set of volume flow rate and power consumption readings 
have been taken for three blowers / experimentally . 
These readings are plugged in the quadratic equation 
and they have been solved for the constants, so this 
procedure finally results m the form 


3.429 (Vol^^ 2 + 1.143 (Vol air > + 3 ?» 67 


‘ blower 


, kW 


1000 x ri 


m 


f ^*18) 


where / 


vol 


air 


r\ 


m 


3 

The volume flow rate of air/ m /s 

The mechanical efficiency of the fan (0,85) 


3.7 EVAPORATIVE COOLING PROCESS 

figure 3.5(a) exhibits the evaporative cooling 
arrangement for the room. The unsaturated outside air 
comes in contact with the wetted surface where the evapo- 
ration of the water lowers the humidified air temperature. 
This process is schematically represented by process AB 
in 3.5 (b). The temperature of the humidified air is 
related to the desert cooler efficiency (^p) as 



T db - (T <Sb ' T «b > 


(3.19) 


T = dry bulb temperature of the outside air/°C 
db 4 

rh -j wet bulb temperature of the outside air y °C 
wb 

T sa dry bulb temperature of the air leaving 
the desert cooler# °C • 


where/ 
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The volume of humidified air for the airconditioning is 
calculated from* 


d,air 


5 TR h"“T 

air a,o a,i 


P 

air 


1.026 


(T 4 - 


T- ) 

dg 


where P is the density of the humid air. 

a! Jl 


( 3.20) 


3.8 CONSTRAINTS FOR THE PROBLEM* 


As mentioned earlier/ rour variables, being restricted 
between some practical values, are taken into consideration 
for the optimization. So, these constraints are mathematically 


expressed 

as 

follows: 



4.0 

<_ 

X(D < 

15.0 

3.21( a) 

2.0 

< 

X( 2) < 

5.0 

3.21(b) 

1.0 

< 

X( 3) < 

4.0 

3.21(c) 

1.0 

< 

X( 4) < 

4.0 

3.21(d) 


The constraints are transformed into equality constraints 
by taking the lower and upper limits, separately. As 
a result, they are transformed in the following form. 


GG( 1) 

= 4.0-X(l) 

3.22( a) 

GG( 2) 

= X(l) - 15.0 

3.22( b) 

GG( 3) 

= 2.0 - X( 2) 

3.22( c) 

GG( 4) 

= X(2) - 5.0 

3.22(d) 
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GG( 5) 

= 1.0 - X( 3) 

3.22C e) 

GG( 6) 

= X( 3) - 4.0 

3.22(f) 

00(7) 

= 1 o - X( 4 ) 

3.22(g) 

GG( 8) 

= X( 4 ) - 4.0 

3.22(h) 


After normalising the equations/ the constraints are 
considered for the optimization of the problem. 


3.9 ESTII1ATI0N OF INITIAL COST 

As mentioned in Sec. 3.6/ the design cooling load 
is calculated for Kanpur 7 in May as m [9] . After evalu- 
ating these loads, the equipments are selected on the basis 
of their capacities. In this regard, a due consideration 
has been given to purchase the equipments on the basis of 
recommended capacity ranges available in the market. A 
stepwise variation is adopted for the costs and capacities 
of the various equipments. This feature is very well 
incorporated in the computer progranme. 

Thus, the total initial cost for the system is 
expressed as 


Cl - ( c CO mp + ^evap + c cond + ^Tslower^ X 

C 3.23) 


where, 

c s= cost of the compressor, Rs 
comp 

C = cost of the evaporator, Rs 
evap 
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c cond “ cost of the condenser, Rs 

C blower = cost the blower, Rs 

MF = Multiplication factor to take care of 
the miscellaneous cost (1-2) 

L = life of the system, years. 


of 

When a hybrid system ^air conditioning and desert cooler 
is adopted between the months April and June, if feasible, 
the initial cost, Cj, is considered to be higher and 
given by, 


C I C I + Sesert^d ^ 3 * 

where , 

C_ = initial cost of the airconditioning system, 
Rs/year 

r , = cost of the desert cooler, Rs 

^Desert 

Tj. = life of the desert cooler , years, 
d 

In the present work, both these cases are studied 
separately and results are compared. 


3.10 estimation of running cost* 

For the running cost, the electricity cost is the 
primary factor to be considered. Keeping m view, the 
variation m electricity cost, a futuristic statistical 
approach has been adopted to make the whole formulation 
as a generalised algorithm [13 ] . Proceeding as m [13], 
the electric charge, ( C Q ) ,is found to bes 


24) 
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1.2S3 

0.97986 + exp (-0.09338 I) 


( 3.25) 


where, I refers to the duration of the years from the 
present one. 


The effective electricity cost is given by [15] 


°E 


1 

L 


L 

S 

1=1 


(1+ R ) 1 " 1 


( 3.26) 


where / 

C = cost of electricity m the ith year ,Rs/kWh 
R = interest rate 

L = life of the machineries, years 

So, the running cost for compressor (c^) is given by. 


C 


1 


comp 


x Gg, Rs/h 


( 3. 27 ) 


For the evaporator fan 7 


C 2 " R fan x V Rs/h 


For the condenser fan/ 


c 3 = P fan 31 S' Rs/h 


( 3.28) 


( 3.29) 


The running cost of all these accessories is calculated 
on the hourly basis by taking into account by sinusoidal 
variation of outside air temperature and summed upto get 
the running co-st/year. The maintenance cost is taken 
to be 10% of the initial investment [15]. 
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'k 

Mathematically, it can be expressed as , 

c, = (t.l) c T LktEZ ... - 1 — (3.30) 

1 R( 1+R) L - 1 

The total running cost is given by 

+ C 2 + C 3 + c 4 ^ Rs/year (3.31) 


In the case o± combined operation of airconditioning and 
desert cooler the running cost includes for the desert 
cooler only for the feasible period during April to June* 


Then the total running cost be come s. 


** 


C R ~ C E * D,C ratmg * ( ' C 1* C 2 +C 3* + C 4 


( 3,32) 


where C„ the effective electricity cost, Rs/kwh 

D.C ^ The desert cooler rating, kw. 

rating 

Thus, the total cost of the system, Qp,is given by 


Qp = C R + C x , Rs/year 


(3.33), 


In Eq.(3.30) when desert cooler is also used, C-j. is replaced 
by Cj- of Eg. (3.24) . 

* Thin operative cost is for the period during which, it 
operates alone. 
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In the case of hybrid system, it is 

I 

C T = C R + C I / Rs/year 


C 3,34) 


For generalisation purpose, the equations 3^33 and 3.34 
are modified as. 


Cnp ( C„ + C_) 

__ - — , kWh/year 

E C B 


3.35C a) 





(C R + C T ) 

-p, , kWh/year 

a 


3.35(b) 


Thus the ratio QpAV, can be written in functional form as, 

Qp/Cg « f(X 1 , X 2 , X 3 , X 4 ), where X ± ' s are the variables 

taken for the present study to get the optimum values for 
the minimum of the ratio Q^/Cg. 

In the present analysis, all these possible combinations 
are incorporated m a generalised computer programme for 
optimization and the same is incorporated m Appendix— G» 
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CHAPTER 4 

RESULTS AND DISCUSS ION 


This chapter deals with the calculation of cooling 
loads on the basis of the theoretical and actual temperature 
variations having the effect of heat capacity and ti selag 
factor of the structures The economic criterion has been 
adopted for the optimum choice of the air conditioning systems 
using R~12 and R-22 refrigerants. A comparison of a hyorid 
system (evaporative cooling and airconditioning systems) 
with a conventional airconditioning system is discussed 
under optimum conditions. For estimation of the optimum 
cost, Qp/Cg, the life of the airconditioner and the desert 
cooler are taken as 20 years and 10 years, respectively# 

The following are the results obtained m the present work. 

4.1 THEORETICAL AND ACTUAL TEMPERATURE VARIATION: 

Figure 4.1 gives the actual temperature variation 
with time as well as the suggested sinusoidal variation in 
temperature based on maximum and minimum temperatures of 
the day. The actual temperature is slightly higher from 
8 A.M. to 3 P M* than the theoretical value before the peak 
hour. Thereafter, the actual temperature is lower than 
that of the theoretical value. The areas under both curves 
reveal small difference ( about 4 to 7%) . For the airconditioning 
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system, the working hours have been taken to be 8 A.M. to 
* 

5PM An underestimation of cooling loads oy using 
sinusoidal variation is found in this. This deviation is 
calculated for three years, 1982, 1983, 1985 and the same 
are tabulated in Tables 4*la to d, for different indoor 
conditions. From these values, the average deviation is 
calculated as 6%. Sc a multiplier 1.06 can be used to 
get the actual total running loads for 10 hours. 

The variation in the sol-air temperatures for 
different directions is shown m Fig. 4.2. It indicates 
that the maximum sol-air temperature occurs on the western 
facing wall at 3 P.M. 

4.2 EFFECT OF HEAT CAPACITY AND TIMELAG FACTOR ON THE 
COOLING LOADS. 

Figures 4.3(a) and 4.3(o) represent the effect of 
the heat capacity and timelag of the structures on the 
hourly cooling load estimation. From this, it is found out 
that the estimation of the cooling load without taking into 
account of heat capacity of the structures leads to an 
overestimation of 20.29% than that of the cooling load 


The computer programme is quite general and can be used 
for any duration. 
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calculated for the structure with the heat capacity. This 
difference is quite pronounce. When the wall thickness 
was halved, the above is found to be 18.02%. 

4 3 OPTIMUM RESULTS FOR VARIOUS INDOOR CONDITIONS i 

For the indoor conditions of 24 °C DB and 65% RH, 
the design cooling load is 9.8% and 24.5% higher than 
that for the indoor conditions of 27«5°C DB, 56% RH 
and 30°C DB, 60% Rh, respectively. Figures 4.4 and 4.5 
represent the variations m initial , running and total 
costs with respect to the rise in condensing temperature 
(T^; of the system, for R-22 and R-22 refrigerants, 
respectively. In this case, the approach temperature is 
taken as 10°C, while the other design variables, the fall 
m lower temperature At^ the degree of superheat A T sup , 
and the degree of subcooling AT gub are maintained at 
the optimum values. The optimum total cost^/Cg, for 
R-12 refrigerant system is 4.752% less than that of R-22 
refrigerant system for the same operating variables* 

Figure 4.6 shows that the optimum total cost, 
(Qp/Cg), for a hybrid system of evaporative cooling and R-12 
air conditioning system yields 4% gam over the conventional 
air conditioner used from April to October, every year* This 
is applicable for the indoor conduction of 27.5°C dry-bulb 
temperature and 56% RH, with approach temperature as 10° C. 
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Figure 4.7 depiqrts the comparison of the optimum 
total cost between a conventional airconditioning system 
and the system with evaporative cooling. The optimum total 
cost tor the conventional airconditioning system is 5.562% 
higher than that the system used with evaporative cooling 
facility. The variation m the optimum total cost (Cp/Cjg) 
with respect to the fall in the lower temperature of the 
system, At^ is displayed in Fig. 4. 8. The effect of AT^ 
on the optimum total cost is less than that of A T^« 

Figures 4.9(a) and (b) snow that the degrees of superheat 
and subcooling have relatively little effect on the optimum 
total cost of the system. 

Figure 4.10 represents the variation m the 
optimum total cost with respect to the approach temperature 
of the systems for various indoor conditions. The approach 
temperature is almost proportional to the optimum total 
cost or vice versa. The various values of optimum variables 
for R-12 and R-22 refrigerant systems are given in the 
Table 4*2, along with the optimum total cost, penalty function 
values and the penalty parameters. 



TABLE * 4.2 . OPTIMUM QUANTITIES FOR R-12 AND R-22 AIRCONDITIONING SYSTEMS 


>i 

M I U 
rfl rtf (D^* 

<3 rtf 0 ) U 

& Q) S'-- 


C 
>i O 
+J H 
H +> <tf 
tf O 3 
c; ch 
5 ) 3 <0 
Qi IH > 


e o h 

•H H P 
4 J 4 J \ 

aji b 

O M u 


27 37 $7 £7 37 
« M S w S w H « n H 


m S *? 1:0 «? tHCOCO fO in ro 
Sri S A H * fclHl CD I 


■B 

n O 
£h o 


£ 

co O 

fa ° 

< 


rH U 

&h o 
< 


XJ 

Eh O 
<3 ° 


CM ^ 

ON CD 

rH \D 

r-l ON 

r* r* 

• » 

o o 


rH ON 

cm co 

O ro 

vD O 


CO CM 

ON ON 

O o 

o o 


CM 

LO 

co 

CM 

ro 

m 

ON 

ro 

ro 

rd 

co 

ro 

d» 

ro 

CM 

d 1 


r- 

ON 

00 

# 

• 

• 

• 

o 

O 

o 

O 

CM 

o 

CO 

m 

rH 

r* 

i — 1 

CM 

• 

• 

■ 

• 

tH 

o 

d 4 

co 

O 

rH 

ro 

vO 

00 

LO 

CM 

I> 

o 

r- 

CO 

r- 

CM 

00 

in 

O 

CM 


CM 

in 

A 

• 

• 

» 

tH 

O 

d* 

CD 

o 

rH 

CO 

D 

CO 

m 

CM 

r* 



CD 

r- 

m 

rH 

ro 

rH 

ON 

o 

co 

rd 

CM 

ON 

d* 

ON 

d* 

ON 

d* 

r* 

f 

* 

• 

• 

CM 

rH 

CM 

rd 

H 

d< 

d* 

tH 

CM 

in 

m 

rH 

O 

rH 

id 

d* 

O 

o 

o 

id 

• 

4 

* 

ft 

rH 

CM 

*d 

rd 

m 

rH 

CM 

LO 

CD 

ON 

vO 


d* 

ON 

d* 

rH 

rH 

CO 

O 

O 

• 

• 

• 

* 

CM 

tH 

CM 

CM 

rH 

o 

CO 

o 

CD 

CM 

€0 

rH 

r- 

r- 

CO 

d 1 

in 

in 

CM 

rH 

• 

• 

• 

ft 

in 

d< 

d 4 

d 1 


rtf u 

o a; a) 

M a M u 

04 & 30 

o, (U +> 

4J <0 


d 4 vO 


X) d< 
rH 


o o 

tH tH 


I 

u 

a) 

ON 

d 

*H 'O 
MH -U 4 ) 
0) C W 
04 03 O 


s s s s is 

.S •5 *5 • E • B 

+> O +> O -P U +>U 4J U 

H • H • d * d • H • 

so so £ P & p so 


+J 4-> +J 
3 3 3 

B .B .S 


4J 

° § 
si jc 


+J O +> U -P U -P -P U -P 

H H d d d d 

SOSPSPS 3:0 & 


CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

tH 

rH 

rH 

rH 

CM 

rH 

rH 

rH 

rH 

CM 

CM 

I 

I 

I 

1 

1 

1 

1 

1 

1 

1 

1 

« 


04 

P 4 

05 

05 

05 

05 

04 

05 

04 


O W 

a S 
LO 

• ^ 


CONTINUED 


60 




cdCt/Ce i kWh/Yrco 









10 3 l x10 3 l xio 3 


For indoor conditions of 27.5°C DB and 56%RH 

Refrigerant R-12 
^approach g 10°C 

SYSTEM WITHOUT DESERTCOOIER 

SYSTEM WITH DESERTCOOLER 


7.5L 3.8L 3.2L 


0 4 8 12 16 

AT h , °C 
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Fig. 4.10 Variation in optimum total cost with approach 
temperature. 
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CHAPTER 5 


CONCLUSIONS AMD SUGGESTIONS 


5.1 CONCLUSIONS: 

i) A multiplication factor has been found out for 
the estimation of actual cooling load with 
the help of maximum and minimum temperatures 
of a given location. 

ii) The heat capacity and time lag factor render the 
cooling load about 20% less than that of 

without these factors. 

iii) The evaporative cooling system can be substituted 
for the conventional air-conditioner during the 
period of April to June in the case of inside 
design condition of T db = 30°c, 0 = 60% and 

T db = 27 5 °C, 0 = 56%. 

iv) The implementation of the hybrid evaporative 
cooling and airconditioning system is found to 
be more economical than the conventional 
airconditioner, under optimum conditions. 

v) For the inside conditions of 27.5°C D3, 56% RH 
and 30 9 C D3, 60% RH, the suggested hybrid 
system yields a gain of 4 % and 5.6% over the 
conventional airconditioning system, 
respectively • 
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vi ) The opt i-aum total cost of the comfort 

airconditioning using R— 12 is cheaper than 
that of R— 22 system by 4.75% for the inside 
design conditions of 24 °c D3 and 65% RH. 

vii) The variations in condensing and evaporator 

temperatures have more influence on the total 
cost of the system than the degrees of sub- 
cooling and superheating. 

viii) The variation in total cost of the airconditioning 
system is almost proportional to the approach 
temperature. 

5*2 SUGGESTIONS: 

i) a combination of evaporative cooling and conven- 
tional airconditioner can be fabricated as a 
single unit. The overall cost can be estimated 
on the basis of actual cost for commercial 
venture* 

ii) The cost of the building may also be included 

in the optimization in order to make the approach 
even more general* 

iii) This approach can also be extended to the 
vapour-absorption system under optimum conditions* 
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APPEND IX -B 


Table -3 

VALUES FOR SOLAR INTENSITY CALCULATIONS 


— 

t 

Month, ! 

i 

i 

i 

i 

i 

Equation of 
t i me/mi nu t e s 

~i 

• 

| Declination 

1 of angle 

t degrees 

« - . 

j A 1 

! W/m 2 

1 3 1 

t 

I 

1 

I 

c i 

Jan. 

-11.2 

-20.0 

1230 

0.142 

0.058 

Feb. 

-13.9 

-10.8 

1214 

0.144 

0.060 

March 

- 7.5 

0.0 

1185 

0.156 

0.071 

April 

1.1 

11.6 

1135 

0.180 

0.097 

May 

3.3 

20.0 

1103 

0.196 

0.121 

J une 

-1.4 

23.45 

1088 

0.205 

0.134 

July 

—6 « 2 

20.6 

1085 

0.207 

0.136 

August 

-2.4 

12.3 

1107 

0.201 

0.122 

September 

7.5 

0.0 

1151 

0.177 

0.092 

October 

15.4 

1 

h- 4 

o 

» 

c* 

1192 

0.160 

0.073 

November 

13.8 

-19.8 

1220 

0.149 

0.063 

December 

1.6 

-23.45 

1233 

0.142 

0.057 


= Apparent solar irradiation at air mass = 0 , W/m' 
= Atmospheric extinction coefficient 
C 1 « Diffuse radiation factor. 
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APPENDIX - c 

EQUATIO NS FOR THE BSTIMAT ION OF PRIMARY ANGLES 

The equation for, LST, local solar time 

LST = LCT 4 - Equation of time (C-l) 

where^ 

LCT = 1ST + ( 82*5-Longitude of the 

location) /I 5, hrs (c-2) 

The equations for the hour angle a^, are given by 
a h = (LST - 12.0) .15/ for LST < 12.0 

(C-3) 

= (12.0 - LST). 15, for LST > 12.0 

The relationships between the latitude angle 1, the hour 
angle a. , declination angle d, altitude angle j3 , and 
azimuth angle y are. 

Sin P =( Cos 1. Cos h. Cos d + Sin 1. Sin d) (C— 4) 

Cos y = Sin (Cos 1. Sind - Cos d. Sin 1. Cos (C— 5) 

For wall solar azimuth angle a , 

a = [ ^ -Y I, for east and west facing walls, 

= Y / for north facing walls. 


= f(Y-lt)l, for south facing walls. 


(C-6) 
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The relationship for the incidence angle , © is 

Cose = CcosP) ( Cosa) (Cosj<y + ( s ln , (sln for 
tilted surfaces, with tilt angle jZf^ 


= ( CosP) (cosa), tor verti cal surfaces 

= Sir p,- for horizontal surfaces. 
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APPEND XX ~D 


EQUATIONS FOR THERMAL PROPERTIES OF GLASS 


glass. 


The relationships for the thermal properties of the 

transmittivit y ref lectivity jc ^and absorptivity, a f 

9 g 

t r = ( 1-r 1 ) 2 . a c /(l - A a 2 ) (D-l) 


where , 


( 1 - r* ) 


r = 


= 1 


(1 - 


r ' - 


,2 

r* . a c ) 

/ - . 1 v 2 

(1 - r ) 


(1-r* . a c ) 


r ‘ = 


, S±n (©, - © ) 

1 r 3 S_ 

2 L 


Sin 2 (©_ + © ) 

g g 0 

Sin" © 

“ E cV (1 “ 2 ^ 

exp C g n 2 


tan (© 


tan (©_ 


where 



Angle of incidence 

Angle of refracted rays 
Extinction coefficient 
Thickness of the glass (m) 
Refraction index* 


(D-2) 


(D-3) 



] 
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BUILDING PARTICULARS FOR COOLING LOAD CALCULATIONS 

Room Size = 8.4 x 4.0 x 3.8, m 3 

Thickness of the wall = 0.4 m 

Thickness of the roof = 0.2 m 

Number of windows; 

On the Northern wall = 2 
On the Southern wall = 2 
On the Eastern wall = O 
On the Western wall = O 

2 

Area of the windows = 1.08 m 

Number of doors: 

On Western walls = 1 
On the other sides = O 
Number of occupants = 25 

Number of fans inside the room = 4 
Number of lights inside the room = 8 
Ratings: 

= lOO watts 


Fans 

Lights 


40 watts. 



88 


APPEND IX-F 



89 


APP "2ND IX — F 

PROPERTIES OF REFRIGERANTS 


Refrigerant R-12: 

c p (T) = 0.59524 + 0.00181715 (T+15) 

s f (T) = 0.142093 + 0.33747 (T/100) - 0.0399 606( T/100) 2 

+ 0.019868( T/100) 3 + 0.0140 23 (T/lOO) 4 
hg( T ) = 188.86 + 0, 44027 8(T) - 7 .02007 ( T/lOO) 2 

- 5.07651 (T/100) 3 - 3.82545 (T/100) 4 


h f (T) = 36.1554 + 0.9 2fcl0S(T) + 6.89916 (T/100) 2 

+ 3.73414 (T/lOO) 3 + 5.91673 (T/lOO) 4 


Ref rig rant R-22: 

c p (T) = 0.60 + 0.0005246 (T+4 5+0.04 3686 (T+45) 2 )), 
for -60 < T < 10 

= 0.70114 + 0.0029529 ( T-IO+O. 01047 2( T-IOJ 2 ) 
for lO < T < 50 

S f (T) = 0.181332 + 0.43667 (T/100) - 0.0596756 (T/lOO) 2 

O jt 

- 0.0234343 (T/100) + 0.0638047 (T/100) 

h (T) = 251-06 + 35.4577 (T/lOO) - 19.4993 (T/100) 2 
g 

- 7.62005 (T/100) 3 - 11.67 56 { T/lOO) 4 

h f (T) * 46 - 2102 + 1.20394 (T) + 6.828 (T/lOO) 2 

- 8 -99387 (T/100) 3 + 20,0612 (T/lOO) 4 

where/ T .in °C, c in kj/fog .°C, s f in kJ/kg.°K/ h in kJ/kg. 
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